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I. INTRODUCTION 



Bogomolny-Prasad-Sommerficld monopoles are static solitons occurring in certain (3+l)-dimcnsional gauge field 
theories. They have attracted interest continually since they were discovered over two decades ago. The simplest BPS 
monopoles are SU (2) monopoles, which have an associated integer referred to as the topological charge. Monopoles 
with topological charge k are called fc-monopoles, and their total energy is proportional to k. For any given k there 
are fc-monopole solutions which resemble k well-separated 1-monopoles. 

It is interesting, but often difficult, to examine monopoles associated with gauge groups larger than SU{2). Such 
monopoles have features not found in the SU(2) case, for example, the existence of spherically symmetric multi- 
monopoles 0. 

There is a powerful approach to monopoles; the ADHMN construction. To perform this construction a nonlinear 
differential equation, called the Nahm equation, must be solved and its solution, the Nahm data, used to define a linear 
ordinary differential equation. This linear equation, which we shall refer to as the ADHMN construction equation, 
must then be solved. Its solutions yield the fields via an integration procedure. Recently Platonic symmetries have 
been exploited to construct Nahm data and these symmetric Nahm data have been used to examine some particular 
examples of SU(2) monopoles. In this paper we discuss two cases where the same Platonic Nahm data, slightly 
modified, can be used to study monopoles associated with larger gauge groups. 

Section || is an introduction to SU (IV) Nahm data and monopoles. We show that for minimal symmetry breaking 
the Nahm data for some multi-monopoles is simpler than the corresponding SU (2) Nahm data. This allows us, in a 
simple way, to modify known SU(2) Nahm data so that it is SU(N) Nahm data. It is possible to use SU{2) Nahm 
data to produce SU(2) monopoles embedded in an SU(N) theory. Such embedded monopoles behave like 577(2) 
monopoles. The modification we consider is more radical than a simple embedding and the corresponding monopoles 
behave quite unlike the way SU(2) monopoles do. 

In Section [II we use charge three Nahm data with tetrahedral symmetry to construct a one-parameter family of 
ST/ (4) monopoles with tetrahedral symmetry. The dynamics of slow moving monopoles is approximated by geo desic 
motion in the moduli space of solutions Q . Since the one-parameter family of solutions described in Section III is the 
fixed point set of the action of the tetrahedral group, it must be a geodesic. Thus the one-parameter family described 
in Section III is an example of a pathological scattering process in which the monopole never separates into distinct 
objects. 

In Section IV we use charge four Nahm data with tetrahedral symmetry to construct a two-parameter family of 
ST/(5) monopoles with tetrahedral symmetry. This two-parameter family is totally geodesic in the whole moduli 
space. Under the assumption that the transformation between the metric on the space of Nahm data and the metric 
on the moduli space of monopoles is an isometry, we undertake a numerical study of the low energy dynamics of 
SU(5) tetrahedral monopoles. We find an exotic dynamics involving both single and double scatterings through 
configurations with octahedral symmetry. 



II. MONOPOLES AND NAHM DATA 



BPS SU(N) monopoles are topological solitons in an SU(N) Yang-Mills-Higgs gauge theory with no Higgs 
self-coupling. They are finite energy solutions to the Bogomolny equation 

A$ = -\e ijk F jk (2.1) 

where D; = + [Ai, is the covariant derivative with Ai an su(IV)- valued gauge potential and Fjk the gauge field. 
$ is an su (IV)- valued scalar field, called the Higgs field. Non-trivial asymptotic conditions are imposed on the Higgs 
field, which are responsible for the existence of topological soliton solutions to the theory. It is required that, as 
r = |x| approaches infinity, $ takes values in the gauge orbit of the matrix 

M = idiag(^i,^ 2 , • ■ ■ ,Hn)- (2.2) 

By convention it is assumed that \x\ < /12 < ■ ■ • < /ijv- Since <& is traceless fi% + /i2 + . . . + = 0. This M is the 
vacuum expectation value for $ and the symmetry group of M under gauge transformation is called the residual, or 
unbroken, symmetry group. Thus, for example, if all the \i v are different the residual symmetry group is the maximal 
torus [/(l)^ -1 . This is known as maximal, or generic, symmetry breaking. The soliton solutions are associated with 
IV — 1 integers; this is because the boundary condition on <I> implies a map, $00, from the large sphere at infinity into 
the quotient group 
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$oo : S 2 -> orhiW (w) M = SU(N)/U(1) 



N-l 



(2.3) 



and 



7T 2 (SU(N)/U{1 



,w-n 



7Tl(C/(l 



,JV— 1 



JV-l 



(2.4) 



In contrast, this paper concerns the minimal symmetry breaking case, in which all but one of the n P are identical, so 
the residual symmetry group is U(N — 1). It is convenient to choose /ii = — (N — 1) and fi p — 1 for p — 2 . . . N. Since 



7T2 (SU (N) /U (N — 1)) = Z 



(2.5) 



there is only one topological integer associated with solutions. Nonetheless, a given solution has N — 1 integers 
associated with it, these arise in the following way; careful analysis of the boundary conditions indicates that there is 
a choice of gauge such that the Higgs field for large r, in a given direction, is given by 



$(r) = zdiag(^i,/i 2 , . . . ~ 7T di ag(&i, k 2 , . . .,k N ) + 0(r 2 ). 

2r 

In the maximal symmetry breaking case the topological charges are given by 



(2.6) 



(2.7) 



In the case of minimal symmetry breaking only the first of these numbers, mi, is a topological charge. Nonetheless, 
the remaining m p constitute an integer characterization of a solution. This characterization is gauge invariant up 
to reordering of the integers k p . The m p are known as magnetic weights, with the matrix diag (k\, k 2 , ■ ■ ■ , few) often 
called the charge matrix. 

There are some obvious ways of embedding su(2) in su(N), for example, 



a (3 
-0 -a 



a ... (3 



-(3 ... -a 



V 



(2.8) 



Important SU(N) monopoles can be produced by embedding the SU(2) 1-monopolc fields, which are known su(2)- 
valued fields, in su(N). Some care must be taken in producing these embedded monopoles to ensure that the 
asymptotic behaviour is correct; the SU(2) monopole may nee d to be scaled and it may be necessary to add a 
constant diagonal field beyon d the plain embedding described by (2.8), details may be found in Obviously there 



is an embedding of the form (2.8) for each choice of two columns in the target matrix. The embedded 1-monopoles 
have a single k p = 1 and another k p = —1, the rest are zero. The choice of columns for the embedding dictates which 
two k p are non-zero. Recall that in the case of minimal symmetry breaking the choice of order of the k p is a gauge 
choice. In fact, in the case of minimal symmetry breaking, the embedded 1-monopole is unique up to position and 
gauge transformation. Solutions with k\=k have k times the energy of this basic solution and so it is reasonable to 
call these k p = k monopoles /c-monopoles. There are of course different types of such fc-monopoles corresponding to 
different magnetic weights. 

Consider SU(3) monopoles with minimal symmetry breaking. For k = 2 there are two distinct types of monopoles, 
those with TO2 = and those with W2 = 1. The mi = 2 case is equivalent to the mi = case; m,2 can be changed 
from to 2 by reordering k 2 and k%. If m 2 = the monopoles are all embeddings of su(2) 2-monopoles and this 
case is not interesting as an example of sit (3) 2-monopoles. The m 2 = 1 case has been studied by Dancer S and by 
Dancer and Leese Q , by considering Nahm data. 

There is an equivalence between Nahm data and BPS monopoles. In the case of SU(N) monopoles the Nahm data 
are triplets of anti-hcrmitian matrix functions (Ti,T 2 ,T3) of s over the intervals (fi p , fi p +i). The size of the matrices 
depends on the corresponding values of m p ; the matrices (Ti, T 2 , T$) are m p x m p matrices in the interval (fi p , [i p+ i). 
They are required to be non-singular in each region and to satisfy the Nahm equation 



dTi 1 r , 

~te = r ijk[Tj > Tk] - 



(2.9) 
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There are complicated boundary conditions prescribed between the Nahm data in abutting intervals, which are 
detailed in Nahm's original paper 0. We will follow Hurtubise and Murray's formulation of the Nahm data boundary 
conditions for distinct \x p |l and then take the limit of coincident \x p to describe the minimal symmetry breaking case. 

For ease of notation we shall describe the case where m p _i > m p (ie. k p < for p > 1) since a similar result holds 
after a reordering if this is not satisfied. 

Monopoles are constructed from their corresponding Nahm data by first solving a first order differential equation in 
which the Nahm data appear as coefficients. This is called the ADHMN construction equation. This choice affects the 
order of the k p . Rather than describe it in full generality, it will be described below in the particular form required. 
The matching and boundary conditions on Nahm data are designed to ensure that the ADHMN construction equation 
has the correct number of solutions required to yield the correct type of monopole fields. Define the function 



JY 



k(s) = k p 9(s - fip) 



(2.10) 



P =i 



where 9(s) is the usual Heaviside function. In the interval ([i p , k(s) = i 
shape depends on the charge matrix of the corresponding monopole. If k(s) near fi p is 



It is a rectilinear skyline whose 



k n 



k(s) 



then as s approaches fi p from below we require 

Ti(s) = 



iRi + 0(l) <3( 2 (IM-i)/2) 

oV fcpl_1)/2 ) t; + o(z) 



where z = s — /i p and where 



as s approaches \i p from above. 



T 4 ( S ) = T[ + 0{z) 



(2.11) 



(2.12) 



It follows from the Nahm equation (2.G) that the \k p \ x \k p \ residue matrices (Ri, i?2, R3) in ( [2.11 ) form a represen- 



tation of su(2). The boundary conditions require that this representation is the unique irreducible |/c p |-dimensional 
representation of su(2). 

In summary, at the boundary between two abutting intervals, if the Nahm matrices are to p _i x m p _i on the left and 
m p x m p on the right an m p x m p block continues through the boundary and there is an (m p _i — m p ) x (rn p _i — m p ) 
block simple pole whose residues form an irreducible representation of su(2). The boundary conditions for to p _i = m p 
are given in, for example, |8]. While these boundary conditions are involved, their function is simply one of limiting 
the number of solutions to the ADHMN construction equation. 

The 1-dimensional representations of su(2) are trivial. Thus, if k p — —1 for all p > 1, fc(s) is 



V>2 



~ 1 1 
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and the Nahm data has only one pole, it is at s = fX\. It seems reasonable to suppose that this result holds in the 
limit of coincident Thus, if we fix \i\ — — (N — 1) and [i v ^\ = 1, we expect that the N x N Nahm data whose sole 
pole is at s — — (N— 1), satisfying the Nahm equations and having acceptable residues, are the Nahm data of SU(N) 
monopoles with minimal symmetry breaking. For this to be the case it is only required that the ADHMN construction 
equation have the correct number of solutions over the interval. This index calculation is easily performed using the 
methods of ||. The topological charge of the corresponding monopole solution is, of necessity, k\ = N — 1, since the 
k p must add to zero. The magnetic weights are each one less than the proceeding one. We say that the magnetic 
weights are distinct. It has recently been proven by Nakajima that all monopoles of this type can be constructed from 
the described Nahm data 0). 

The Nahm data for su{2) fc-monopoles are k x k anti-hermitian matrix solutions of the Nahm equation. They have 
poles at s — — 1 and s = 1. It is obvious that this Nahm data, 



k 








k 




i : 





can be used to generate Nahm data for SU(k + 1) fc-monopoles with distinct weights 

1 



k 



1 



In examples where the charge k SU(2) data is known, the SU(k + 1) data is generated by a translation and rescaling 
of s so that a pole occurs at s = fix but the second pole is moved outside the interval s € [jttivMjv] i e - it is lost from 
the Nahm data. The 2-monopole Nahm data is known exactly, and was used by Dancer in || to construct SU(3) 
monopoles. This is the simplest application of the above procedure. Platonic symmetry groups have previously been 
used to derive higher charge Nahm data, and in this paper we discuss the corresponding SU(k + 1) monopoles. 



III. SU(4) MONOPOLES WITH TETRAHEDRAL SYMMETRY 

In the previous Section we discussed the Nahm data for SU{k + 1) monopoles with minimal symmetry breaking, 
charge k and distinct magnetic weights. For the remainder of the paper it will be convenient to perform a translation 
s i ► s , so that a pole in the Nahm data always occurs at s = 0. In this Section we describe some aspects of the 
the ADHMN construction, which calculates the monopole fields ($, A{) from the Nahm data. We then go on to apply 
this construction to obtain a one-parameter family of monopoles with k = 3, which have tetrahedral symmetry. 

Given Nahm data (T\ , T 2 , T 3 ) for a k- monopole we must solve the ADHMN construction equation, for s € [0, k + 1] , 

(l2fe^ + Ife O Xjcr-j + iTj ® cTj-)v = (3.1) 

for the complex 2k- vector v(s), where 1/c denotes the k x k identity matrix, <7j are the Pauli matrices and x = 
(xi, X2, X3) is the point in space at which the monopole fields are to be calculated. Introducing the inner product 

rk+l 

(vi,v 2 > = / vlv 2 ds (3.2) 
Jo 



5 



then the solutions of (3.1) which we requir e are those which are normalizable with respect to ( |3.2| ). It can be shown 
that the space of normalizable solutions to (3.1) has (complex) dimension k + 1. If Vi, .., v^+i is an orthonormal basis 
for this space then the ijth matrix element, of the Higgs field is given by 



(3.3) 



A similar expression exists for the gauge potential, but the energy density, £ , may be computed without calculating 
the gauge potential by using the formula 



£ = Atr($ 2 ) 



(3.4) 



where A denotes the Laplacian on ]R 3 . 

For most of the example s co nsidered in this paper the Nahm data is sufficiently complicated that the matrix 
linear differential equation (3.1) can not be solved analytically in closed form. In these cases we use a numerical 
implementation of the ADHMN construction which involves solving the ordinary differential equations using a standard 
fourth order Runge-Kutta method. The numerical implementation is similar to that introduced by the authors in 
fiof , but is simplified by the fact that the Nahm data we consider here has a pole at only one end of the s interval. 
This eliminates the need for the shooting part of the numerical algorithm described in [ flO[ . 

In references pl] , |l0| , ^2| it is explained how Platonic symmetry (that is, tetrahedral, octahedral or icosahedral 
symmetry) may be applied to Nahm data for SU(2) monopoles of charge k. We use these, as explained in the 
previous Section, to easily obtain the solutions to Nahm's equation for SU(k + 1) monopoles. 

From pOl] we have that the Nahm data for tetrahedrally symmetric monopoles with k = 3 has the form 



(3.5) 
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where 



2p(ws) p{ws) 



and p is the Weierstrass function satisfying 



P~ 



where ' denotes differentiation with respect to the argument. 

For SU(2) monopoles the boundary condition requires that the Nahm data has a simple pole at s 
determines n to be 



± K — Kq 



r(i/6)r(i/3) 

12v^ 



(3.6) 



(3.7) 



2 which 



(3.8) 



For SU (4) monopoles the corresponding boundary condition is less restrictive, namely that we require the Nahm data 
to have no poles for s £ (0,4]. Given the SU(2) result this implies that k must satisfy the condition 



kq/2 < k < Kq/2. 



(3.9) 



It is a simple task to show that the remaining requirements for Nahm data are satisfied, and hence we have proved 
the existence of a one-parameter family of SU (4) monopoles with tetrahedral symmetry. The one-parameter family 
is given by n in the above interval and the corresponding family of spectral curves is 



r/ 3 + i36K 3 C(C 4 - 1) = 0. 



(3.10) 



Note that k = gives the spectral curve rf = 0, which is the spectral curve of a 3-monopole with spherical symmetry. 
Such a spherically symmetric monopole was constructed several years ago W by using a spherically symmetric ansatz 
in the field equations. We shall now see how this solution arises in the ADHMN construction by explicitly calculating 
the Higgs field in this case. 

Taking the limit k — > and using the property that 



Pit- 



as u 







(3-11) 



G 



gives 



z = -1/s 



(3.12) 



in this limit. 

Since the monopole with this Nahm data is spherically sym metri c, we only need to constructing the Higgs field 
along an axis. Setting (x±, x%, X3) — (0,0, r) the linear equation (3.1) becomes 



rfv 

ds 



(r-l/s) 










-(r-l/s) 
-y/2/a 











-y/2/8 
r 












-r -y/2/s 

-y/2/a (r + l/s) 

















-(r + l/s) 



(3.13) 



Clearly the first and last components of v decouple and are elementary to solve. The remaining four equations may 
be decoupled into two pairs, which can then be converted into two second order equations and solved by a Laplace 
transform. The full regular solution is 

cise~ rs 
ca^pl-e"" (l + £)) 



c 3 (- 



sinh(r 



e rs (l-^)) 



-c 3 V2(— ^ - 



C4.se 



(3.14) 



where c\, 02,03, C4 are arbitrary constants. If we select the orthono rma l basis where v-i has only Ci non-zero, then 
the Higgs field will be diagonal. Performing the required integrals (3.3) gives the result 



($) 



33 



. 6r - 3 + e" 8r (64r 3 + 48r 2 + 18r + 3) 
~ % 2r(l-e- 8j -(32r 2 + 8r + l)) 
e 8r (-4r + 3) - 384r 3 + 64r 2 - 40r - 6 



(3.15) 



~(128r 3 + 128r 2 + 44r + 3) 



•(-4r+ 1) + 128r 3 - 8r - 2 



,-Sr 



(128r 3 + 64r 2 + 12r+l) 



with ($)22 and ( < i ) )44 obtained by the replacement r i~+ — r in (&)n and ( ( &)33 respectively. It is a simple task to 
verify that indeed this solution has the correct asymptotic behaviour ie. 



$ -> idiag(-3, 1,1,1) 



(3.16) 



To compute the Higgs field for non-zero values of k is a much more difficult task, since all the components of the 
vector v become coupled together. Thus we turn to the numerical implementation of the ADHMN construction. 
In Figure 1 we display the results in the form of three dimensional plots of surfaces of constant energy density for 
the values k/k = —0.25, —0.10, 0.00, +0.10, +0.25. As the parameter k increases from zero, the spherical monopole 
deforms into a tetrahedral monopole. As k — * ±«o/2, the monopole approaches the embedded SU(2) tetrahedral 
3- monopole asymptotically. In fact, even for the value k — —kq/A (Figure 1.1) the energy density looks very similar 
to that of the SU(2) tetrahedral 3-monopole jl(|. Note that changing the sign of k gives a monopole corresponding 
to the dual tetrahedron. 
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FIG. 1. Tetrahedral scattering of an SU(4) 3-monopole. 

In the moduli space approximation P| the dynamics of k monopoles is approximated by geodesic motion on the 
fc-monopole moduli space Mk- From the spectral curve approach it is clear that, having fixed an orientation and 
centre of mass, we have constructed the unique one-parameter family of tetrahedrally symmetric SU(A) 3-monopolcs. 
Since the fixed point set of a group action gives a totally geodesic submanifold, this one-parameter family is a geodesic 
in M 3 . Hence, within the moduli space approximation, this family of monopoles may be interpreted as describing the 
low energy dynamics of three deforming monopoles. Using this interpretation we see from Figure 1 that during the 
course of the motion a tetrahedron gets smoothed out into a sphere which then deforms back into the dual tetrahedron. 
This gives an example of dynamics in which the monopoles never become separated. 

It is known in the case of SU(2) monopoles that there are closed geodesies In the geodesies approximation 

such geodesies correspond to periodic monopole motions during which the monopoles never separate. That is not the 
case here, here the motion is not periodic and the geodesic is not closed; it runs between points on the asymptotic 
boundary of the moduli space which do not correspond to separated monopoles. 

Obviously the method applied in this Section to tetrahedral SU(4) 3- monopoles can easily be carried over to 
construct SU(k + 1) /c-monopoles, given the Nahm data for an SU(2) fc-monopole. In general, given a p-dimensional 
family of SU(2) monopoles there will be a corresponding (p + l)-dimensional family of SU(k + 1) monopoles. Thus, 
for example, it is a simple task to construct the Nahm data for the one-parameter family of octahedrally symmetric 
SU(5) 4-monopoles which derive from the unique octahedrally symmetric SU(2) 4-monopole p~l] ] . However, it is more 
interesting to consider geodesic motion in the two-dimensional moduli space of SU(k + 1) monopoles derived from a 
one-parameter family of ST/ (2) monopoles corresponding to a geodesic in the SU(2) moduli space. Physically, this 
will allow us to examine how the dynamics of SU (k + 1) monopoles compares with the dynamics of SU (2) monopoles. 
We shall do this in the following Section, for the case of SU (5) 4-monopoles with tetrahedral symmetry. 
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IV. SU(5) MONOPOLES WITH TETRAHEDRAL SYMMETRY 



After fixing the orientation and centre of mass, there is a one-parameter family of tctrahedrally symmetric charge 
four SU{2) monopoles @. The associated Nahm data takes the form 



Ti(s) =x(s)X i +y(s)Y i 
where the tctrahedrally symmetric Nahm triplets are 



z(s)Zi 



( 



(X 1 +iX 2 ,X 3 



(Y 1 + iY 2 ,Y 3 ) = 4 



(Z 1 +iZ 2 ,Z 3 ) = V3 2 
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The reduced equations for the three real functions x, y, z can be solved to yield 



/ \ K ( /—, — r 1 p'( K s) 



1 p'(ks) 



z(s) 



2p(ns) 



Here p is the Weierstrass function satisfying 

p' 2 = 4p 3 - 4p + 12a 2 

with prime denoting differentiation with respect to the argument. 

The spectral curve for tetrahedrally symmetric 4-monopoles has the form 

rf + im ? C(C 4 - 1) + /? 2 (C 8 + 14C 4 + 1) = 



(4.3) 

(4.4) 
(4.5) 

(4.6) 
(4.7) 



where a and [3 are real constants. The relation between these constants and those appearing in the above Nahm data 
is given by 



36a«T 



(3 2 = 3k\ 



(4.8) 



In the SU(2) case, the requirement that the Nahm data has a second pole at s = 2 means that n must be taken to 
be half the real period of the elliptic function (4.6). Thus, k is determined given the parameter a, and we have the 
required one-parameter family. Furthermore, a is restricted to lie in the interval a G (— a c , a c ), with a c = 3- 5 / 4 V2. 
The elliptic function becomes rational at a — ±a c , with infinite real period so that there is no second pole, and hence 
there is no corresponding SU (2) monopole. 

Applying the boundary conditions for SU(5) monopoles is a different story: we now require no singularities of the 
Nahm data for s £ (0, 5]. If we consider a £ (— a c , a c ) then the result is similar to that of the previous Section. The 
r ang e of k is now restricted, k £ (— kq, ko) where 5ko is the real period of the elliptic function (4.6). Using the formula 
(4.8) this determines a domain in the (a, (3) plane of the spectral curve coefficients. For a = ±a c , there is no second 
pole of the elliptic function, so the value of k is unrestricted. This case corresponds to two curves in the (a, (3) plane, 
which pass through the origin and continue off to infinity. 

To examine the case \a\ > a c we need to con side r some properties of elliptic functions p6| . For \a\ < a c the 
discriminant of the elliptic curve determined by (4.6) is positive and the period lattice is rectangular. The elliptic 
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function has poles on the real axis, but no zeros. However, for \a\ > a c the character of the elliptic function changes 
since the discriminant is now negative. The period lattice is rhombic an d in add ition to having poles on the real 
axis, the elliptic function also has zeros on the real axis. From equations ( fi~3| - |4.5| ) we see that a zero of the elliptic 
function also corresponds to singular Nahm data. Thus, in this case, there is a restriction on k given by k G (— kq, kq) 
where 5/5o is the smallest real root of the elliptic function (4.6), that is, p(5ko) = 0. This defines a second domain in 
the (a, (3) plane which matches smoothly onto the first, with the joining boundary being the curves determined by 
a = ±a c . 

We now have no restriction on the parameter a, so we must also consider the limit a — > oo. In this limit it can be 
shown that kq — > 0, but in such a way that the combination okq is finite, though it can be non-zero. In terms of the 
spectral curve constants this limit corresponds to monopoles with ft = 0, but a restricted only to lie in some finite 
range. We refer to such monopoles as purely tetrahedral, since the octahedral term in the spectral curve is absent. 
This is an interesting result, since no such purely tetrahedral charge four monopoles occur in the ST/ (2) theory Of 
course, given the existence of purely tetrahedral monopoles it is simple to study the reduced Nahm equations directly 
in the case (3 — and obtain the same result as above without the need for limit taking. 

From the above analysis it is seen that in order to compute the Nahm data and calculate the domain of definition in 
the (a, pi) plane, numerical algorithms must be employed to compute not only elliptic functions and their derivatives 
but also their periods and elliptic logarithms. In the SU(2) case this task was much easier, since for a rectangular period 
lattice the required computations can be performed using Jacobi elliptic functions with real arguments. However, in 
the rhombic case this is not true, and it is better to work directly with the Weierstrass function. Standard algorithms 
are used which are based upon the AGM method and truncated series [17[ . 
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FIG. 2. Five geodesies for SU(5) tetrahedral monopolcs. 



In Figure 2 we plot (dashed lines) the boundary of the spectral curve coefficients in the (a, (3) plane for a € [—2, 2]. 
Note that we allow (3 to be negative, even though the points (a, (3) and (a, — (3) give the same spectral curve and are 
hence gauge equivalent. The reason is that the coefficient of the octahedral term in the spectral curve is non-negative, 
so that [3 2 is the correct parameterization, rather than, say, (3. This arises since a cube is inversion symmetric, 
whereas the tetrahedron is not and would lead to a change of sign in the spectral curve coefficient. Thus the analogue 
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of the SU(A) tetrahedral geodesic of the last Section, where the tetrahedron deformed to a sphere and out through 
the inverted tetrahedron, is the SU (5) octahedral geodesic along the line a = 0, where a cube deforms into a sphere 
(a,/?) = (0,0) and out through the inverted cube (which is gauge equivalent). 

It should be stressed that our X-shaped representation of the moduli space is not a reflection of any metric properties 
of the moduli space. As mentioned above, we know that the line a = is a geodesic (by symmetry), but in order to 
determine more general geodesies we must first compute the metric on this two dimensional moduli space, and then 
solve the geodesic equations of motion. We shall do this using numerical techniques, though we must also make an 
assumption, as follows. 

For SU (2) monopoles it is known that the transformation between the monopole moduli space metric and the metric 
on Nahm data is an isometry [^8[ . However, for general SU(N) gauge groups it has not yet been proved that this 
transformation is an isometry, although this is widely believed to be true. There is also circumstantial evidence, for 
example, assuming this result leads to monopole metrics which reproduce conjectured metrics based upon asymptotic 
knowledge |lS) ] . To make progress we shall assume that this transformation is an isometry, and compute the metric 
on the Nahm data. This assumption was also made in previous studies on the dynamics of SU(3) monopoles 

The scheme to compute the metric is similar to the 5/7(2) case [|o|, to which we refer the reader for a more detailed 
discussion. Note that recently this SU(2) metric has been computed exactly and in closed form This was then 
used to demonstrate the excellent accuracy of the numerical algorithm |2(J. It is likely that the method of could 
also be used in this case to calculate the SU(5) metric exactly, if required. 

The tangent space is computed by solving the linearized Nahm equation 

^ = e iifc [T i ,V fc ] i = l,2,3 (4.9) 

for the tangent vector (Vi, V 2 , V3) corresponding to the point with Nahm data (Ti, T2, T3). Given two tangent vectors 
Vi, Wi, the metric on Nahm data is 

< V h Wi >= - / Vtr(Wi) da. (4.10) 
Jo i=i 

From the tetrahedral symmetry of the Nahm data it follows that the tangent vectors are tetrahedrally symmetric 
so we may write 

V i =q 1 X i + q 2 Y i + q 3 Z i i = 1,2,3 (4.11) 
where q = (qi, q 2 , 33)* is an analytic real 3-vector function of s € [0, 5]. In terms of q equation (|4.9|) is 



M q where M 



Ax -96y -122/5 
— 6y — 16y— 6x — 6z/5 
-Az -32z -Ax - 32y 



(4.12) 



This ordinary differential equation has a regular-singular point at s — 0. Analysis of the initial value problem at s = 
reveals that there is a two-dimensional family of solutions which are normalizable for s £ [0, 5]. They are given by the 
two-parameter, c — (01,02), family of initial conditions 

q - (0,c 1 s 3 ,c 2 s 2 ) t as s - 0. (4.13) 

Using the asymptotic properties of the Weierstrass function we find that the Nahm data has the behaviour 

P 2 s 3 as 2 n 

y~- , z ~ as s ~ 0. (4.14) 

y 120 72 v ; 

Hence to compute the tangent vector dual to the coordinate a requires the choice c = (0,1/72), whereas to 
compute the tangent vector ^ dual to the coordinate (3 requires c = (/3/60,0). The metric can then be computed as 

d d 
oa oa 

92 = < wh > (415) 

53 = < h I > (4 - 16) 
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with corresponding Lagrangian 



,da. , dB , 9 . da . , dB . 



The metric is computed numerically by solv ing e quation (4.12) using a fixed-step fourth-order Runge-Kutta method, 



with the integrations required in e quati on ( 4.10 ) calculated via a composite Simpsons rule. The geodesic equations 
which follow from the Lagrangian ( |4.17 ) are solved using a variable-step Runge-Kutta method, with the derivatives 
of the metric approximated by finite differences. The accuracy of our scheme was such that energy was conserved to 
four significant figures for all computed geodesic trajectories. 

Note from the above that the metric components g 2 and g$ both vanish for j3 = 0, which is a reflection of our choice 
of B 2 as the spectral curve coefficient and implies that all geodesies which cross the a-axis are parallel to the /3-axis 
at the point of crossing. Thus from the numerical point of view we work with the coordinate B 2 when computing 
geodesies, since it is better behaved than B. 

In Figure 2 we show five geodesies (solid lines), three of which pass through the point (a,B) = (1.0,0.4) and the 
remaining two pass through the point (a, B) — (0.0, 0.1). Many other geodesies were also computed, but the qualitative 
features are captured by those shown. Basically, the results show two kinds of scattering that take place. The first 
kind is similar to the SU(2) scattering and occurs when the geodesic does not stray too far away from the SU(2) 
embedding boundary. The four monopoles approach from infinity on the vertices of a large contracting tetrahedron, 
scatter through a cubic monopole, that is, cross the /3-axis, and emerge on the vertices of an expanding tetrahedron 
dual to the incoming one. We show one geodesic of this kind in the upper half plane. The second kind of scattering is 
more exotic and involves a double scattering through a cubic monopole. The remaining four geodesies are all of this 
kind, with three associated with monopoles which approach from infinity with a positive and one with a negative. 
In each case the geodesic first crosses the /3-axis (a cubic scattering) and then crosses the a-axis, instantaneously 
forming a purely tetrahedral monopole, after which it recrosses the /3-axis (the second cubic scattering) and goes off 
to infinity gauge equivalent to the incoming configuration. 

The two types of scattering described above were the only ones found; no geodesies were found with, for example, 
no cubic scatterings or more than two cubic scatterings. In fact the results in this case are similar in spirit to those 
seen in the study of SU(3) 2-monopole dynamics where it was found that up to two 90° scatterings could take 
place. It would therefore seem that this phenomenon of multiple scatterings is the generic situation for general SU(N) 
monopoles with minimal symmetry breaking. 
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